A molecular superfluid: non-classical rotations in doped para-hydrogen clusters 
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Clusters of para- hydrogen (pH2) have been predicted to exhibit superfluid behavior, but direct ob- 
servation of this phenomenon has been elusive. Combining experiments and theoretical simulations, 
we have determined the size evolution of the superfluid response of pH2 clusters doped with carbon 
dioxide (C02)- Reduction of the effective inertia is observed when the dopant is surrounded by the 
pH2 solvent. This marks the onset of molecular superfluidity in pH2. The fractional occupation of 
solvation rings around CO2 correlates with enhanced superfluid response for certain cluster sizes. 

PACS numbers: 36.40.-c, 36.40.Ei, 36.40.Mr, 67.25.dw 



Superfluidity has been well characterized in the 
bulk liquid phase [l]. In a breakthrough experiment, 
Grebenev, Toennies, and Vilesov observed superfluidity 
in helium nanodroplets [2 , introducing a new way of in- 
vestigating finite quantum systems [3 . Only atomic he- 
lium has been known to exhibit this behaviour at liquid- 
like densities, and the direct observation of superfluid- 
ity in a molecular system has remained elusive. The 
most likely candidate for a molecular superfluid is para- 
hydrogen (PH2), because of its bosonic character, low 
mass, and weak inter molecular forces. In order to con- 
firm superfluidity in molecular hydrogen, an experimen- 
tal determination of the so-called superfluid response or 
non-classical rotational inertia [1 is required. 

The spectroscopic observation of nearly free molecu- 
lar rotation in helium nanodroplets [U [3] has motivated 
several studies aimed at providing an understanding of 
superfluidity at the nanoscale in finite systems. Rovibra- 
tional spectra of linear molecules such as OCS in helium 
droplets are characterized by narrow lines indicative of 
coherent rotation and decoupling from the solvent 0. 
Most importantly, the spectroscopic rotational constant 
5 |4] of the molecular dopant is renormalized by the he- 
lium environment. B is inversely proportional to the ef- 
fective moment of inertia of the system, so in a normal 
classical system, the inertia should grow monotonically 
with system size N. A key experiment on OCS-doped 
clusters with a few {N = 1 — 8) helium atoms has shown 
that their effective inertia could actually be greater than 
that of larger droplets [5 . Thus a decoupling mechanism 
must exist in order to stop the growth of the effective in- 
ertia as the cluster grows beyond = 8. The turnaround 
marks the onset of superfluidity [5j. 

Theoretical analysis has later shown that the size evo- 



lution of B was related to that of the superfluid frac- 
tion, defined as the deviation from a classical response 
to rotation |6tt8j . B constants were also used to estimate 
the 'experimental' superfluid fraction [7[ |8]. Theorists 
have predicted superfluidity in clusters [9], but the 
direct observation of superfluid response in doped pll2 
clusters has remained elusive. Simulation work predicted 
superfluidity for OCS-doped clusters [IQl [11] with 
a turnaround of 5 at A/" = 14. However, this predic- 
tion is not yet confirmed by experiment, since published 
measurements only extend to A^ = 7 [12 . In other exper- 
iments, superfluidity was inferred for doped pll2 clusters 
embedded in helium nanodroplets based on the behav- 
ior of the Q branch of the rovibrational spectrum [[T3l - 
[T6] . Moore and Miller performed related experiments on 
doped HD clusters fTTHTQ]- These results motivated our 
search for a system where the superfluidity of pH2 could 
be discerned. 

Here we report confirmation of molecular superfluid 
behavior in doped pH2 clusters, based on spectroscopic 
determination of their non-classical rotational inertia 
probed by CO2. Our main result confirms earlier the- 
oretical predictions that pR2 clusters could become su- 
perfluid |9j, and lends weight to experiments that have 
shown evidence of superfluidity in doped pll2 clusters 
embedded in helium nanodroplets p^HTB] . We have mea- 
sured infrared spectra and performed finite-temperature 
path-integral Monte Carlo (PIMC) simulations to sup- 
port our findings. 

The experimental apparatus has been described else- 
where [20^. Backing pressures in the range of 5 to 25 at- 
mospheres were used, and the jet nozzle was cooled (-20 
to 100 ^C). Effective cluster rotational temperatures were 
0.3 to 1 K. Since the measurements were performed for 
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i^C^^Os and i^C^^Os, the vibrational frequency shifts 
and B rotational parameters reported here were scaled 
to values appropriate to ^^C^^02 for comparison with 
theory. Due to CO2 nuclear spin statistics and the low 
temperature of the jet expansion, only three transitions 
are generally observed for each cluster: i^(0), i^(2), and 
P(2). We fit these lines in terms of three parameters: 
band origin, B rotational constant, D centrifugal distor- 
tion constant. 

Details of the PIMC calculation method are described 
elsewhere [71 [2TH23 . It relies on the so-called worm algo- 
rithm [241 ES] to account for bosonic exchange. We used 
a recent H2-H2 interaction potential [2ff along with our 
own H2-CO2 potential energy surfaces The calcula- 
tions were performed for a low temperature (0.5 K), with 
256 translational and 128 rotational time slices. A total 
number of 10^-10^ Monte Carlo steps for each simula- 
tion yielded satisfactory error bars. Figure la is a mon- 
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FIG. 1. (color online) (a) Observed spectrum with as- 
signed R{0) transitions labeled by cluster size, N. The three 
strongest (off-scale) transitions are the P(4), P(2), and R{0) 
transitions of the CO2 monomer, (b) Calculated chemical po- 
tential (in wavenumbers; 1.0 cm~^ = 1.438769 K) versus N. 
(c) Measured (filled circles) and theoretical (open triangles) 
shifts of the CO2 asymmetric stretch vibrational frequency 
versus N. 



tage of 4 traces where the experimental conditions were 



optimized for each cluster size range. For this reason, 
and because laser power varies, relative line intensities 
are not quantitative. However, their qualitative trend is 
real, reflecting the behavior of the chemical potential (the 
energy difference between adjacent cluster sizes) curve in 
Fig. lb, with N = 5 and 12 being relatively stable (larger 
negative chemical potential) and = 6, 9 and 10 less so. 
Line assignments for < 8 are straightforward, and con- 
firm the accuracy of the theoretical vibrational shifts. For 
N > S the line assignments are not so obvious. However, 
by using the theoretical CO2 vibrational shifts in Fig. 
1(c) to help guide the assignments, a convincing scheme 
is achieved which is more complete and self-consistent 
than any other. The special nature of = 12 is directly 
evident without any reference to theory: its transitions 
are noticeably stronger, indicating a local maximum in 
the cluster population, which is consistent with the local 
maximum in the absolute value of the chemical potential 
in Fig. 1(b). 

Magic numbers indicative of particularly stable struc- 
tures are observed for = 5, 8, 12 14, and 17. In Fig. 
1(c), the subtle change in slope at A/" = 5 is associated 
with filling of a 'donut ring' of 5 pH2 molecules around 
the CO2 axis. For 5 < A' < 9 the additional pR2 start 
populating a second ring. The excellent agreement be- 
tween the theoretical and experimental vibrational shifts 
attests to the quality of the interaction potentials used 
in the simulations. 

The structural evolution is shown in Fig. 2 by the 
patterns of the pH2 density contours in the frame of the 
CO2 molecule. Occupation numbers of 'donut rings' and 
'terminal caps' are indicated. The occupation numbers, 
which add up to A/", were obtained by integrating the pH2 
density in each ring and rounding to the nearest simple 
fraction. The boundary region between two rings was set 
as the minimum in the density. A central 'donut ring' is 
filled for N = 5 and a full solvation shell is observed for 
N =17] the latter cluster size marks a more pronounced 
change in slope in the vibrational shift curve of Fig. Ic. 
The filling of a second centred ring starts at A^ = 18. 

Figure 3a shows the effective B constant, 5eff as a 
function A^. Calculations based on a dipole-dipole corre- 
lation function can be compared directly to experiment. 
In earlier work on helium clusters doped with CO2 [22j, 
N2O [7112T], and HCCCN [8 , the calculated B was in 
very good agreement with experiment. The present the- 
oretical values agree remarkably well with experiment, 
and capture the overall behavior of the size evolution of 
^eff^- However, some deviations are observed for N > 13. 
These discrepancies may be attributed to the finite (non- 
zero) temperature of the simulations. 

A key question is whether or not B^^^ really probes 
superfluidity. For systems undergoing rotation, the su- 
perfluid fraction is defined in terms of the non-classical 
rotational inertia using a two-fluid model [28l [29] which 
assumes that the total density of the fluid is the sum of 
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FIG. 2. Contours of two dimensional projections of the pH2 
density in the CO2 frame for representative N values. Whole 
and fractional occupation numbers are shown. 



contributions from normal and superfiuid components: 
p = Ps -\- Pn- The superfiuid fraction is = ^ while 
the normal fraction is = defined as fn = yy, the 
ratio of the effective (/eff) over the classical (/ci) inertia 
of pR2- The calculated fs shown in Fig. 3b corresponds 
to the perpendicular response with respect to the CO2 
axis. The experimental fs was obtained using the ap- 
proach proposed in Ref. [7 , according to which B^^^ 
defines the total moment of inertia of the system via the 
relation B^^^ = ^^exp . The effective inertia of is 

total 

obtained by subtracting the moment of inertia of CO2, 
^eff^ = ^totai ~ ^C02 ' experimental estimate of fs is 
presented in Fig. 3b. The experimental and theoretical 
fs values agree fairly well in their overall behavior over 
a large range of sizes. 

We conclude that within linear response theory, CO2 
can be used as a probe of fs in pH2 clusters. It is also 
possible to use the computed fs to obtain two-fluid B 
constants. Such two-fluid values are presented in Fig. 
3a. Their behavior is very close to that of B^^^^ and 
further confirms that CO2 is a probe of superfiuid re- 
sponse. To highlight the non-classical behaviour of /eff, 
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FIG. 3. (a) Beff versus N from experiment (filled circles), and 
dipole-dipole (open squares), two- fluid (open triangles), and 
classical (open circles) calculations, (b) Experimental (filled 
circles) and calculated (open triangles) fs versus N. 



the effective B constants based on classical cluster iner- 
tia {Bel = 2(7 i-f/cQ ) ) ^^^^ shown in Fig. 3a. As 
expected, B^ decreases monotonically with N. 

Figure 4 shows the low and high pH2 density regions in 
3D for two cluster sizes [30 . For = 12, the maximum 
in /s, the low-density surface is continuous, except at 
the ends of the CO2 axis. This is consistent with the 
enhanced bosonic exchanges associated with its larger fs . 
For TV = 17, the low-density ring regions are disconnected 
from one another, consistent with a lower fs. For A" = 12, 
the high-density surface is still highly delocalized within 
each ring, although the rings are no longer connected. 
For A" = 17, pronounced localization is observed at high 
density, and one can clearly enumerate the individual p}l2 
molecules. The orthographic top view for A" = 12 (Fig. 
4b) shows that the central ring is swollen to accommodate 
an additional fraction of a particle, while a similar view 
for A" = 17 (Fig. 4e) shows that all three rings have 
similar radii with staggered particle placement. 

Using CO2 to probe superfluidity for doped pB.2 clus- 
ters in the A^ = 1 — 18 range, we provide the first direct 
measurement of fs in a molecular system. A maximum 
in fs is observed at A" = 12. This value is confirmed 
by consistent experimental and theoretical results, and 
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FIG. 4. (color online) Three-dimensional pH2 densities [30] , 
(a) Side view for N = 12. (b) Top orthographic (or flat) 
view for N — 12. (c) Perspective view for N — 12. (d) side 
view for N = 17. (e) Top orthographic view for N = 17. 
(f) Perspective view for N = 17. The green (dark) and grey 
(light) colors respectively represent high and low densities. 



is associated with the appearance of fractional ring occu- 
pation. Beyond this size, fg decreases due to enhanced lo- 
calization. There is a fine balance between the forces that 
cause localization, and bosonic exchanges that favour su- 
perfluididty. Localization dominates at higher N. A 
full solvation shell with pronounced localization is ob- 
served at A/" = 17; the computed fs there is substantially 
lower than the experimental value. A possible explana- 
tion for this discrepancy is the finite (non-zero) temper- 
ature character of the PIMC simulation and the fact this 
fs estimate is based on linear response theory. A lower 
temperature simulation could lead to a higher fs due to 
quantum melting j3TH33] . Whether or not larger clusters 
will be superfluid remains an open question. The lo- 
calization modulated superfluid response may be viewed 
as the nanoscale analogue of the insulating to superfluid 
transitions observed in cold gases j34] . 
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